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1. Introduction: ME vs. PS 

The production of multi jet final states is one of the searching grounds for new physics 
at current and future collider experiments. Therefore their Monte Carlo simulation, both 
for signal and background processes is of crucial importance for the success of the experi- 
ments. Currently, two alternative ways to fill the phase space for particle emission can be 
formulated: 

1. One might use exact matrix elements (ME) at some given perturbative order in the 
coupling constant (s), say a^. In moxst cases, these MEs are available at the tree level, 
i.e. at the lowest perturbative order for the production of a specific final state. Then, 
the final state particles are usually identified with jets and their phase space is cut 
accordingly. In some cases higher perturbative orders are known as well, and blur this 
simple correspondence. However, in this paper we will focus on tree level MEs, there- 
fore we can stick to the simple identification of partons with jets. The virtue of the 
MEs is that they are exact and take all interference effects at the perturbative order 
into account. The downside is that more and more diagrams have to be considered 
with rising numbers of outgoing particles populating a multidimensional phase space. 
Clearly, our calculational abilities do not live up to this situation. Furthermore, for 
the transition of the partons into hadrons a phcnomenological fragmentation model 
has to be invoked with parameters, that arc a priori free and depend on the fragmen- 
tation scale. Hence, applying a fragmentation scheme immediately at the hard scale 
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of the ME necessitates a re-tuning of its parameters for each cm. energy in order 
to gain reliable results on the hadron level. This is an impossible task and limits the 
applicability of using pure MEs considerably. 

2. Alternatively one might use the parton shower (PS). By expanding around the soft 
and collinear limits of parton emission the radiation pattern factorises into individual 
parton branchings, each giving rise to potentially large soft and collinear logarithms. 
Via multiple emission according to the individual branching probabilities the PS then 
resums all the large logarithms. In fact, the coherent, angular ordered PS is correct 
up to Next -to Leading Logarithmic (NLL) accuracy, i.e. it correctly resums all terms 
of the form a" (L'^" + L-^""^) , where L denotes a large logarithm of the form In s/Qg 
with s the hard scale of the process and Q'q the soft parton resolution scale. Obviously, 
the PS is able to connect both the hard scale and the fragmentation scale. Hence, 
using the PS allows us to tune the fragmentation parameters at one scale, say at 
LEPl energies, where huge statistics were amassed, and to use the same parameters 
at all other scales. However, due to its factorised structure the PS misses important 
interference effects that may play a significant role for detailled analyses. 

To gain some insight into the effect of these interference effects, let us compare the emission 
of a gluon in e^e^ — > qq. Pictorially we can write: 



ME over PS 




Figure 1: Pictorial representation of contributions to single gluon emission in e+e^ qqg as 
given by the ME and the PS. As indicated, the PS does not take into account the interference 
contributions present in the ME. This leads to the ratio ME /PS as depicted in the contour plot. 
At the soft and collinear boundaries of the phase space in the X1-X2 plane (where Xi corresponds 
to the energy fraction of the quark and antiquark) the PS correctly reproduces the ME, whereas in 
the region of hard gluon emission the PS omits the (destructive) interference contribution. 

From the figure above, it can be seen clearly, that in the ratio dayiF,/ da-ps the omitted 
interference terms have a large impact in the region where the gluon is hard and emitted 
at a large angle. Obviously it would be of great interest to combine both ME and PS 
to take full advantage of their respective virtues. A number of attempts in this direction 
were made in [2]- [11], some of them also incorporate NLO corrections for specific processes. 
For the latter, a recent paper [12] provides a general method to construct process specific 
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algorithms to combine MEs at NLO with the PS. This method, however, is limited by the 
number of processes that are calculated at NLO and that the user wants to implement. 
In contrast, in the framework of this paper I would like to restrict the discussion on the 
merging of MEs at tree level with the PS. This allows for a process independent method 
that can be applied immediately. In fact, what I propose in this paper is the extension of 
a previous paper [1], addressing the merging of arbitrary MEs for jet production at tree 
level and the PS in e^e~ annihilations, to processes with hadronic initial states. 



exact ME 
L0 4jet 



NLL 

resummed 
in PS 



^act ME 
LO 5jet, but also 
NLO 4jet 



The basic idea of this method is to divide 
the phase space for parton radiation in two 
disjoint regions, one region of jet production 
that is filled with help of MEs and a com- 
plementary region of jet production, that is 
filled with the PS. The combined ME+PS 
then should correctly reproduce all large log- 
arithms up to NLL accuracy and the sub- 
leading terms present in the ME. Consider 
as an example the process e+e" — > jets. The 
respective orders in and L are given in the 
figure to the right. For the case of four jet 
production, the method incorporates all the 
contributions labelled with a 4, i.e. the full 
tree level ME, plus all NLL contributions of 
higher perturbative order, represented by the 
green blobs to the right of the vertical four 
jet line. 

The challenge in such a procedure is to avoid 
double counting. In our example of four jet 
production this translates into ensuring that 
the green blobs labelled by a 4 are taken into 
account only once. This can be achieved by 

applying weights on the tree level ME such that the NLL terms of higher perturbative order 
in Og are correctly treated and by a veto on the production of extra jets in the PS. To see 
how this works, let us stick first to the example of jet production in e^e~ annihilations 
before going to more general cases. 




Schematic picture of perturbative orders 
in e~^e~ — >jets. Clearly, for every order 
in as one additional gluon can be emit- 
ted that can be soft or collinear leading to 
maximal two more large logaritms. 



2. ME-I-PS in e+e annihilations 
2.1 Jet rates at NLL accuracy 

I'd like to begin the discussion of the combination procedure with a brief review of the 
PS. There, individual branchings, i.e. parton emissions, are governed by the Sudakov form 
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factor, 



A(r, t) = exp < 



T z+it') 

It! 



(2.1) 



yielding a probability for no branch of parton a into partons h and c between scales 
T and t. In Eq.(2.1) the scales T and t might denote virtual masses like for instance in 
PYTHIA [13] or scaled opening angles like in HERWIG, [14]. z denotes the energy fraction 
of the parton after branching, the limits on the ^-integration depend on the meaning of 
the scale parameters as does the particular form of the transverse momentum in the 
argument of ag, and the splitting function Pa-^hc depends on the parton types involved. In 
the following we will consider only angular ordered, coherent PS [15]. Then the scales T 
and t play the role of scaled opening angles, i.e. the branching scale of parton a, ta reads 



and the z-limits are given by 



^t^/ta <Z<1- ^to/ta 



(2.2) 



(2.3) 



where to is connected to a minimal branching angle. 

However, up to NLL accuracy we can then perform the ^;-integrals in Eq.(2.1) yielding the 
NLL-Sudakov form factor 



A^^^(Q,g)=exp 



with the integrated splitting functions 



dq'Vq,g{q',Q) 



(2.4) 



2Cf as{q) A Q 3 

rg(g,Q) = log — - 

TT g \ g 4 

2CAas{q) / Q 11 
r,(,,Q) = --^^log--- 



(2.5) 



In Eqs.(2.4) and (2.5) I have explicitly denoted the parton type quark (q) or ghion (g) - 
and the scales Q and q are now directly related to transverse momenta. For further details 
I would like to refer to [16]. 

Now we are in the position to give expressions for n jet rates in e^e~ annihilations in the 
Durham- or k±_ scheme [17, 18]. In this scheme two partons belong to different jets, if 



2min{£;2,£;2|(i _ ^^^q^ 



(2.6) 



where the parameter j/jet regulates the "hardness" of the jets and s = is the cm. 
energy squared of the e'^e~ pair. Defining the jet resolution scale 



Q. 



jet 



(2.7) 
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we have 



NLL 



NLL 





(2.8) 



Qjet 



The interpretation is quite simple: Since on the parton level jet rates are mere probabilities, 
the two jet rate is just the combined probability that neither the quark nor the anti quark 
have emitted another parton at scales above the jet resolution scale. In turn the two jet 
rate is a combination of two possible "histories". In either history, one quark does not emit 
a parton resolvable at Qjct, whereas the other one first propagates down to an intermediate 
scale q without having radiated a parton resolvable at Qjet- Then, at q it decays into a quark 
and a gluon and both decay products experience an evolution down to the jet resolution 
scale without branching any further. That way, the rates for higher jet configurations can 
be constructed as well, for details see [18] . However, let us note that omitting the integral 
over q differential jet rates at NLL can be constructed, i.e. expressions like dTZs/dq. 

2.2 Differential jet rates and the ME weight 

This leads us directly to the construction of the weights on the ME. Having chosen the 
final state momenta pi according to the ME with ag taken at the jet resolution scale the 
weight is construct along the following lines: 

1. Construct the correct "PS history": 

• Merge the two partons i and j with the smallest yij. Combine only "allowed" 
pairs, like for instance {quark, gluon} or {quark, anti-quark}. The momentum 
of the new parton is the sum of the momenta of i and j. 

• Repeat the step above until only a qq pair remains. 

2. For each parton line of type p between gin and gout apply a weight 



4. Accept or reject the kinematical configuration according to the combined weight. 




(2.9) 



3. For each QCD node apply a correction factor 

OisiQjet) 



(2.10) 
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As an example let us consider three jet production as in 
Fig. 2. 2. In this case the correction weight reads 




^q{Q}et,q)^g{Qiet,q) 



yielding the correct jet rate TZs after taking into account 
configurations where the gluon gets clustered with the 
other quark line and after inserting the approbriate in- 



tegrated splitting function and integration over all scales q. 

In this way, all the integrated splitting functions of the differential jet rates at NLL accuracy 
have been replaced by the correct ME. This leads to an improved description of the region 
Uqg , yqg > yjet, i-©- the region of hard jet emission where interference effects matter more 
and more. Clearly the combined weight of this procedure and the ME takes into account 
the full ME for jet production at tree-level plus all the higher order leading and next-to 
leading logaritmic contributions. 

From the considerations above we can formulate the re-weighting procedure as a cookbook 
recipe as follows: 

1. Expand the PS weight up to the perturbative order in for jet production at tree 
level. In fact this boils down to a product of integrated splitting functions. 

2. Replace this leading term of the result with the full ME squared. 

Let me note in passing that this is exactly the approach utilised by Pythia [13] and Herwig 
[14] but only up to the first order in q^. In both programmes, either the first or the hardest 
emission off both the quark and the anti-quark line is re- weighted with the ME ^, thus 
suitably replacing the splitting function. The re-weighting procedure makes use of the 
fact, that dau^iqqg) < daps{qqg) in the full phase space of gluon emission. This relation 
is not valid for higher parton configurations any more and consequently the re-weighting 
procedure has to be modified. The most obvious way to keep the re-weighting applicable 
is to enhance the value of ag in the Sudakov form factor, ag ^ k ■ ag with A; > 1 a suitable 
constant factor. However, due to a potentially large number of rejected emissions, this 
modification might lead to tremendously inefficient parton showers. For more details on 
the rejection procedure in various processes see [5]- [7]. 

2.3 Vetoing the PS 

Having corrected the ME on differential jet rates exact up to NLL accuracy, we are now 
left with the task of performing the PS. The question then naturally arises, at which scale 
to start the PS. Naively one might try to start the parton shower at the jet resolution scale, 
i.e. at Qjet- However, this is wrong. A naive counter argument is that in such a case, there 

^the weight is given by daME{qqg) / dapsiqqg) 



-6- 



would be a radiation dip at scales just below Qjet- For a more formal argument, consider 
two jet events. Their rate is 

(2.11) 

see Eq.(2.8). The probability for a two jet event at an even smaller resolution scale Qq < 
Qjet should read 

T^2''HQo) = [Af L(£:^^,g^j)]2 . (2.12) 

If we start the PS at the jet resolution scale, however, the probability for two jet events at 
scale Qo would read 



which is clearly wrong. The correct answer is to start the shower at 
the hard scale Ecm and veto all emissions with q > Qjet- Identifying 
vetoed emissions with crosses in the figure on the right we then find 
for one quark line 



^q{Ecm, <5o)Ag(£^cm, Qjet) 



1+ J dqTg{q,E^^) + ... 

Qjet 



Ag(-Kcm, Qo)Ag(.Ecm, Qjet) 
A5(£'cm, Qjet) 



Ag(£;cm,Qo),(2.14) 



exactly as demanded. Similar reasoning holds true for more jets, again I want to refer to 
[1]. This leads to the following algorithm for the PS: 

• The PS evolution for any outgoing parton starts at the nodal scale q > Qjet where 
the parton was produced. 

• In any branch at some nodal value q, the softer of the two outgoing partons is pro- 
duced at q whereas the harder parton is produced at some larger scale. This is of 
specific interest for g ^ gg nodes. 

Taken everything together it can be shown that in the combination of re-weighting the 

MEs and vetoing the PS correct NLL jet rates are reproduced at all resolution scales above 
the fragmentation scale. In particular, the dependence on Qjct cancels at NLL accuracy. 
Furthermore, colour configurations are chosen in a gauge independent way. 



3. ME+PS for hadron collisions 

For the construction of a combination procedure with similar virtues wc employ the same 
strategy of dividing the phase space for parton emission into a region of jet production 
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and a region of jet evolution. Again, this division is achieved by means of a A;_l -measure 
[19, 20], and again MEs together with a suitable weight are responsible for the production 
of jets whereas the PS together with a veto on the emergence of extra jets takes care of 
the jet evolution down to the fragmentation scale. 



3.1 Some brief reminders 

Before I go into more detail I'd like to review briefly the PS in the initial state. For 
definiteness let us constrain ourselves on the case of processes with hadrons as initial 

states, since the extension to processes involving both hadrons and leptons in the initial 
state (like for instance Deep-Inelastic Scattering, DIS) is straightforward. First, some 
process /ii/i2— > X and its configuration in phase space is generated according to 

an,h,^x =J2 [ dxidx2d^f^lixi,Q^)f^^^{x2,Q')^^^^^^ , (3.1) 

where Xi are the energy fractions the partons pi carry. Furthermore, da/d^ is the partial 
differential cross section of the partonic subprocess w.r.t. the differential phase space ele- 
ment $ of the outgoing particles. The scale entering the parton distribution functions 
(PDF) depends on the phase space point under consideration. In Drell-Yan processes, for 
instance, this scale is given by the invariant mass of the lepton pair, whereas in QCD-type 
subprocesses this scale is given by the of the outgoing partons. 

Starting from this scale the PS in the initial scale now proceeds backwards down to some 
infrared scale tQ. This backward evolution is due to the fact that in difference to the final 
state PS here the parton configurations both at the hard and the soft scales are already 
known. The idea behind the backward evolution is to include the PDFs in order to generate 
only physically meaningful parton ensembles during the PS evolution. 
In this framework, the probability for a parton b to be evolved backwards from {t2,X2) to 
{ti,X2) with no branching resolvable at the scale to reads 

''^*'''''"^^-A(x„i.)A,(t„io)- ^'-'^ 

The inclusion of the PDFs actually ensures the correct physical behaviour at low and high 
values of X2- Having chosen a scale ti, the corresponding value z = X2/X1 of the splitting 
a ^ be has to be selected according to 

Ois\PL{z,ti)] X2/zfa{x2/z,ti) „ . x 

27r X2fb{x2,h) 

After each branching, the system is boosted into the cm. frame of the two outermost 
incoming partons. 

Let us continue the discussion now with the outline of the algorithm I'd like to propose 
before I start constructing the weights on the MEs. 
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3.2 Proposed algorithm 

The algoritm I propose goes as follows: 

1. Cluster initial and final state particles backwards according to the longitudinal in- 
variant kj_ scheme [19, 20] until a 2 ^ 2 process remains. This clustering is achieved 
step by step in the cm. system of the incoming partons. In this scheme the initial 
and final state partons are included in the following fashion: 

• If the two particles considered are both outgoing, their measure yij is given by 

2mm{Ef,En{l- cos 9ij) 

y., = . 

s 

see Eq. (2.6). Here, s is the invariant mass squared of the outgoing particles, 
are their transverse momenta w.r.t. the beam axis and 77 and (j) are their 
pseudorapidities and azimuthal angles. 

If two outgoing, i.e. final state, particles are clustered, the resulting particle 
again is a final state particle with p = Pi + pj- 

• If one of the two particles, say j is one of the two incoming partons then 

2Efil-cose.,)^pl^_ (3.4) 

s s 

If an incoming and an outgoing particle are clustered, the new particle is incom- 
ing, and its momentum is p = pj — pi. Note that in such a case a boost is in 
order to the cm. frame of the new pair of two incoming particles. 

In this scheme, the outgoing and incoming particles are clustered "towards" the hard 
2^2 process, i.e. clustering two particles leads to another particle resolved at a 
higher scale. 

Note that in case we consider DIS-like processes the measures are given by the 
energies and the cosines, whereas in case of purely hadronic initial states the measures 
are given in terms of transverse momenta. 

2. Find the hardest in the "core" 2 2-subprocess. Examples for this hardest k 
are: 

• s = Mfi in Drell-Yan type qq — >■ II subprocesses. 

• in QCD subprocesses. 

3. Apply a weight constructed by comparing histories. Some examples will be given 
below. 

4. Start the initial and final state parton showers from the 2 2-subprocess with 
corresponding starting conditions, i.e. every leg starts its evolution at the scale where 
it was produced. 
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5. Veto on jet emissions in the PS. 
3.3 Constructing the ME weight 




As an example let us consider Drell-Yan processes. Similar to the case of e~^e~ annihilations 
we first construct a PS weight, that consists of ratios of Sudakov form factors and splitting 
functions. Taking tjet as the jet resolution scale and t as the hard scale t = Mf^ this weight 
reads for the configuration in the example above 

q{xi/zi,ti) Ag(t,ijet) 



Wp5 



ziq{xi,t) Ag(ti,ijet) 

q{xi/ziZ2,t2) Ag{ti,tjet 



Z2qixi/Zi,ti) Ag{t2,tjet) 
g(xi/ziZ2,ijet) Ag{t2,tjet 



■ as{tl)Pg^qg{Zl) 

■ Ots{t2)Pq^qg\ 



Z2) 



(3.5) 

q{xxlzxZ2M) Aqit-^et.tict)' 
We have to combine this with the PDFs contained in the expression for the cross section 
for pure Drell-Yan processes with Q'^ = t, 



daoY 



dxidx2q{xi,Q^)q{x2, Q'^)da qq^u{s , as{Q^)) 



see Eq. (3.1) and compare the result with a similar expression for the production of two 
additional gluons at a lower scale, 

dUDY+gg = j dx'^dx2q{x'^ = Xi / ZiZ2,^ct)q{x2,t]ct)da qq^iigg{s / Z1Z2, as{^et) ■ 

We can easily read off the correction weight on the matrix element, namely 

as{ti) as{t2) 



Aq{t,^ 



(3.6) 



Ois{t]et) as(tjet) 

Note that the missing factor ziZ2 in the denominator as well as the missing PDF in the 
numerator is compensated by the differential cross section for the "Drell-yan + 2 gluon"- 
process at the jet resolution scale. 
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Of course, things become more complicated, if some flavour changing branchings like for 
instance gluon splitting occurs between the jet resolution and the hard scale. For the 
example above the PS weight would read 



g{xi/ziZ2,t^ct) Ag{t2,tjet) 



(3.7) 



g{xi/ziZ2,t2) Ag(tjet,ijet) ' 

resulting in the following correction weight 

We read off that the correction factor is a combination of Sudakov form factors depending 
only on the nodes of emissions and the flavours of the incoming lines. This is in fact quite 
similar to what we already knew from pure final state considerations. 

3.4 Vetoed show^ers 

Finally, let us turn to the PS. For the final state particles, the PS evolution is done in 
the same fashion as in the case of e+e~ annihilations. For the two initial state lines the 
evolution from the harder - spacelike - scales down to lower scales proceeds as follows: 

1. Evolve to lower scales with help of the Sudakov form factors. 

2. Veto on all emissions with t > tjet- 

3. For the first allowed branching with t^w, < tjet the PDF part of the weight reads: 

>VPDP = 4P;^. (3.9) 

Zf{x,tjet) 

In the two examples above, x = xi/ziZ2- 



4. Summciry and outlook 

In this paper I have suggested a method to combine MEs and PS for hadronic interactions. 
This method is an extension of the already known one for e+e" annihilations [1] that has 
been implemented and successfully tested vs. experimental data in the new event generator 
AMEGIC++/APACIC++ [21, 22]. 

In both cases the method correctly takes into account the full perturbative order in ag of the 
tree-level ME for the corresponding multi jet production. Additionally, all contributions 
up to NLL accuracy are reproduced in the jet rates via a re-weighting procedure. Suitable 
starting conditions and a veto applied on the PS avoid double counting. 
Nevertheless, I would like to stress that the extension of the method I proposed still waits 
for the proof of correctness and for the implementation into an event generator. This is 
work in progress. 
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